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(^-eometrical  Acoustics  I:   The  Theory  of  Weak  Shock  Waves'' 

Jose^oh  B.  Keller 


I .   Introduction 

One  of  the  main  difficulties  in  the  theoretical  analysis 
of  shock  wave  problems  is  that  the  shock  front  and  the  flox-j 
behind  it  must  be  determined  simultaneously  since  they  affect 
each  other.   However,  for  very  weak  (acoustic)  shocks  the 
interaction  is  so  slight  that  the  motion  of  the  shock  front  and 
the  variation  in  its  strength  can  be  determined  indeuendently 
of  the  rest  of  the  flow.   This  fact  is  deduced  in  the  present 
report,  and  the  main  conseauences  of  it  are  also  determined. 
The  result  is  that  weak  shocks  can  be  analyzed  by  the  methods 
of  geometrical  optics,  and  for  this  reason  the  theory  of  weak 
shocks  may  be  called  geometrical  acoustics.   Another  acoustical 
phenomenon- -the  propagation  of  periodic  sound  waves  of  high 
frequency--can  also  be  analyzed  by  the  methods  of  geometrical 
acoustics,  and  this  analysis  will  be  presented  in  another 
report. 

Although  it  has  often  been  noticed  that  certain  shock 
wave  phenomena  could  be  described  in  optical  terms,  there  does 
not  seem  to  be  any  complete  derivation  of  geometrical  acoustics 
in  the  literature.   The  present  derivation  begins  with  a 
deduction  of  the  discontinuity  conditions  for  curved 
discontinuities  in  any  continuous  medliom,  and  continues  with 
the  analysis  of  the  possible  discontinuities  in  arbitrary 
media.   Then  the  acoustic  equations  are  deduced  as  the 
variational  equations  of  the  nonlinear  fluid  dynamic  equations. 
In  the  same  way  the  acoustic  discontinuity  conditions  are 
obtained,  and  it  is  shox^rn  how  the  geometrical  optics, 
considerations  of  wavef rents,  rays,  etc.  apply  to  the 
determination  of  shock  waves  in  perfect  fluids.   Finally  as  one 
of  the  main  results,  the  variation  of  shock  strength  along  a  ray 
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is  determined.   This  rGsn.lt  is  extended  to  the  determination 
of  the  reflection  and  transmission  coefficients  for  shocks 
incident  on  contact  discontinuities  (boundaries). 

For  the  most  effective  use  of  the  present  theory  it  is 
important  to  take  account  of  diffracted  rays  and  wavefronts, 
which  can  bo  done  by  slightly  refining  the  methods  of  ordinary 
geometrical  optics.   An  explanation  of  this  refinement  will  be 
presented  in  a  future  report,  although  use  has  alroadj'"  been 
made  of  it  (see,  for  example,  Reflection  and  Diffraction  of 
Pulses  by  Wedges  and  Corners,  Albert  Blank  and  Joseph  3.  Keller, 
Comm.  Pure  Appl.  Math.,  June  19'^1,  Vol.  IV,  a'o.  2). 

The  present  analysis  is  patterned  after  that  of 
R.  K.  Luneberg's  study  of  the  propagation  of  discontinuities  in 
electromagnetic  theory  (see,  for  example,  his  Mathcaatical 
Theory  of  Optics,  Brown  University  or  his  Propagation  of 
Electromagnetic  Waves,  Now  York  University) . 

1 1 .   Discontinuity  Condi  t  ions  In  Continuuzn  ¥c  chanics 

A.   Fundamental  Equations 

The  fundamental  equations  of  continuum  mechanics,  in 
their  integral  form,  are 

(1)   ^   P^V  +  I  p(u^-v^)dS^  =  0  Mass 

V       S 


2' 


V  V 


'T^,.-pu^(u.-v.  ;  ]dS  .     Linear 
'•'      J   J-     J    Komcntijm 


(3)   ^J  p(e+^u^)dV=j  pf^u^dV  Energy 


V  V 


+  j  [T.,u.-q.-p(e  +  |u2)(u.-v.)]dS. 
S 
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{^) 


^J  o(^i+l"i+2-^l+2^H+l)^^^  Angular 

y  j'lomentum 

V  S 

-  P-^iH-l^l+2-^i+2^i4-inu-v.)]dS. 

(5)  ;^   psdV  -  !  [q.9"^+  ps(u.-v. ) ]dS.  >  0         Second  Law 

V      s 

(6)  e  -   e(o,s)  First  Law 

(7)  e  =  £3 

In  these  equations  S(t)  is  the  surface  of  an  arbitrary 
volume  V(t),  x.  are  Cartesian  coordinates,  i  takes  the  va.lues 
].,2,3  (mod  3)>  p  is  density,  u.  is  Darticle  velocity,  s  is 
entropy  per  unit  volurae,  e  is  internal  energy,  f.  is  external 
force  (both  per  unit  mass),  T..  is  the  stress  tensor,  q.  is 
the  heat  flow  vector  and  9  is  temperature.   The  su'-aiiiation 
convention  is  understood  and  all  quantities  are  functions  of 
X.,  t.   If  the  equation  of  S(t)  is 

(8)  (i)(x^,t)  =  0 
then  the  unit  normal  v.  has  components 

(9)  ^  =  }^^_(^2,)-V2   . 

The  velocity  v.  of  the  surface  is  then  defined  by 

(10)  ^i  = -Mt^'i'x.^'"'^    •      . 

J 

It  is  assumed  that  the  normal  v.  Doints  out  of  V,  and  that 

1  -  ' 

d^.  =  v.dS  where  dS  is  an  element  of  area  of  S. 
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The  preceding  equations  can  be  transformed  into  other 
forms  by  means  of  the  two  transformation  theorems 

(11)  ^J  AdV  =/A^dV  +/Av.dS, 


(12)  |b,   dV^/B.dS-    . 


V   "^i     s 


For  the  validity  of  the  first  theorem  it  is  necessary  that  A 
be  continuous,  and  wo  will  have  to  make  use  of  this  fact. 
Applying  (11)  to  the  left  side  of  (1)  yields 

(13)      J  p^c'^V  +1  pv^dS^  +/  p(u^-v^)d<^^  =  0   . 

V      s        s 

Now  making  use  of  (12)  in  (13)  v;e  find 
(111)  J  [Pt  +  (pui),,j 


V 


,dV  =  0 
1 


Since  ( ll^.)  holds  for  an  arbitrary  volume  V,  wo  have  when  the 
integrand  is  continuous 

(15)  p^  +  ip^i)^    =  0  . 

i 
Proceeding  similarly  with  (2) -(5)  yields 

(16)  u.   +  u.u.   -  p"  T.  .   =  f. 

^t     -'   X  .       ^  ^x .     ^ 
J         J 

(17)  e.  +u.e   +p"-'-[q.  -T.  .u.   ]  =  0 

^    ^  ^1       ^x.   ^J  ^x. 

(18)  T    =  T 

(19)  p(s,  +  u..  s   )  +  (q  9"^),^   >  0 

1  1 


?. 


Equations  (15) -(19)  are  the  differential  forms  of  the 
fundamental  equations  of  continuum  mechanics,  and  apply  only 
when  all  the  functions  p,  u.,  etc.  are  continuous. 

If  any  of  thu  functions  p,  u. ,  etc.  has  a  lump  dis- 
continuity across  a  surface  S' ,  the  above  differential  equa- 
tions do  not  apply  on  S '  .   In  this  case  we  choose  a  rep-ion  V 
which  is  divided  into  two  parts,  V,  and  Vp,  by  s'.   The 
boundary  of  V-,  consists  of  S   and  another  part  which  wc  call  S^  ; 
similarly  Vg  is  bounded  by  s'  and  Sp .   Since  s'  is  a  dis- 
contlnuitjr  surface,  we  indicate  the  liraits  which  functions 
attain  on  approaching  s'  by  a  subscript  one  or  tv70,  according 
as  S   is  approached  through  V-,  or  Vp .   Wc  nov/  apply  (1)  to 
each  of  the  regions  V,  V-j^,  and  Y.^   ^^^^^   then  subtract  the 
equations  for  V-,  and  Vp  from  that  for  V.   ¥c  thus  obtain 


20] 


/  fp(^i-Vi)]idS:i  -  j"  [p(u.-v.)]2dS^2  =   0 


s'  "s' 


Combining  these  integrals,  and  noting  that  ds!-,  =  -ds!p,  we 
obtain 


(21)     /{[p(Ui-v,)]^  -  [o(u,-v,)]2}d 
S' 


ds!o  =  0 


Since  S   i^  an  arbitrary  part  of  the  discontinuity  surface, 
we  have  when  the  integrand  in  (21)  is  continuous, 

(22)  [p(u.-v.)]^v.  =  [p(u.-v.)]2V.   . 
In  exactly  the  same  way  (2),  (3),  (5)  l^ad  to 

(23)  [T,.-pu,(u.-v.)]^v.  .  [T,.-pu.(u.-v.)]pV. 

{Zk)         [T^^.u^.q_.-p(e+  I  u^)(u.-v^.)];l^'j 

=  [1\  ai^-qj.-p(e+  ^  u^  )  (  u^-v  . )  ]2V ^ 
(25)    [q^e>"  +  ps(u^-v^)  ]^V_^  >  [q^©~^+  ps(u^-v^)  jgV^   . 
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In  deriving  (2,9)  it  is  assi-uiied  that  v.  points  toward  side  one; 
if  not  the  inequality  raust  be  reversed.   The  equation  obtained 
in  the  above  manner  from  (k)  is  an  iiMincdiate  consequence  of 
[21].).      Thus  the  discontinuity  conditions  implied  by  (l)-(5) 
are  (22)-(25).   These  conditions  together  with  the  differential 
relations  (15)-(19)  hIso  ii'n/oly  the  integral  equations  (l)-(5). 
(The  derivation  of  (25)  requires  a  proof  that  the  left  side  of 
(5)  computed  for  V  is  greater  than  the  sum  of  the  same 
expressions  computed  for  V-,  and  Vp .   This  follows  directly 
from  (19)  . ) 

3.   Shocks,  Contact  Discontinuities  and  rhasc  Chanf/O  Fronts 
The  jump  conditions  (22) -(25)  may  be  rewritten  in  the 
form 

(26)  p-j^TJ^  =  P2IT2    or   92^^  A  u^  =  -U^^P 

(27)  Pi^iAu.  =  v^AT^^. 

(23)       p^IJ^A(e+|  ut;)  -  v^Z:i(T^.u^-  q.)  =0 

(29)  p-^U-^As  +  v^  A^q^O"-"-)  >  0 

Here  \-jo   have  introduced  the  normal  velocities  TJ-,  ,  Up  of  the 
surface  relative  to  the  medium  on  sides  one  and  tx^ro 
respectively,  and  the  ''jump"  ^,  by 

(30)  ^g    ■=  g-^  -  g2 

(31)  Ui  =  (v.-u,^)v. 

(32^  U  =    (v.-u,,)v. 

2    ^  1   i2 '  1 

If  U,  =0  the  surface  moves  i/ith  the  medium  and  the 
discontinuity  is  called  a  conHta^t  discontinuity.   In  this 
case,  since  p..  and  pp  are  not  zero,  (26) -(29)  become 

(33)  ^2  ^   ^        °^   v^Au^  =  0 

O^l)  V  .  at.  .  =  0 


;  7. 

(35)  ^j^''\1  "  "j-ij^'^i 

(36)  V.  A(q,G"^)  >  G   . 

Equations  (33)  >  (3^)  shov;  that  the  normal  velocity  and  the 
stress  on  the  surfa'ce  are  continuous,  while  (35)  shovjs  that  the 
jurap  in  normal  heat  flu:c  equals  the  jump  in  pov/er  of  the  normal 
stress . 

If  U-,  7^  0,  (27)  may  be  written 


(37)  Au^  =  (p-lU^)"^v^.^T. 


If  in  addition.  A,  p   f^    '.:    the  discontinuity  is  called  a  shock. 
From  (37),  (26)  we  find  U    -i-^g  ^^■ 
relative  to  the  medium  on  side  one 


From  (37)5  (26)  we  find  U    -^^jq   normal  velocity  of  the  shock 


(38)  U,  -  1  1 

"1 


Prom  (38)  >  (26)  we  find  that  Up  is  given  by  (3S)  with  sixb- 
scripts  one  and  two  interchanged. 

If  U,  7^  C  and  A  p  =  0  the  surface  is  called  a  phase 
change  front.   Equation  (37)  applies,  while  (26),  (28)  yield 

(39)  ^i-^^-j  =  '"   I 

v^(u  T.  -q.) 
(1;0)  Tj^  =  -^ -J^^      • 

Equation  (39)  sho^^fs  that  the  normal  velocity  is  continuous, 
while  (I4.O)  gives  the  norntial  velocity  TJ-,  of  the  surface  relative 
to  the  fluid,  which  is  also  eoual  to  Up.   Such  fronts  occur  in 
melting,  freezing,  evaporation  and  condensation,  etc. 

C.   Perfect  Fluids 

A  perfect  (or  non-viscous)  fluid  is  defined  by 

ikl)  T.  .  =  -p5     . 

Here  p  is  the  thormodynamic  pressure,  defined  by 


8. 

(^2)  p  =  -  -^V  =  P^^n    • 

T^or  incompressiblG  media  ( ii.2 )  is  meaningless,  and  then  p  is 
introduced  as  a  new  dependent  fvuiction. 

For  a  perfect  fluid  (l3)  is  automaticallj''  satisfied 
while  (16),  (17) J  (19)  can  be  reduced  to 

(li.3)  u.  +U.U,   +  p'  p   =  f. 

1  ■    :i  1^,    '   ^  X.     1 
ox.         1 

1 


X. 

1 


(■'4-5)  q.  ^V   <  0 

1 

In  deducing  (1-1-5)  j  the  assumption  that  ©  >  0  is  made  use  of. 
At  a  contact  discontinuity  (3l(-)-(36)  become 

(1-I-6)  /\v>   =   0 

ikl)  V  Aq .  =  0 

(I4-8)  q^v.  A(©"^)  >  0 

Thus  p  and  the  normal  heat  flux  are  continuous,  while  heat 
flov7s  from  higher  to  lower  torapcrature  at  a  constant 
discontinu.ity. 

At  a  shock  (3?)j  (3^)  become 

('-1-9)  Au^  =  -Vi(p-L-i)'^Ap 

;"p"p'ap 

(50)  IT  =  i  -£__  . 

At  a  phase  change  front  (k9)  applies,  and  this  with  (39) 
implies  that  u.  and  p  are  continu.ous.   Then  (ii-0)  becomes 

-V  .  Aq  . 

(51)  u  =  ._J._.....I   . 

-■-   p,  Ae 
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If  no  heat  conduction  occnor's,  then  a  phase  change  front  cannot 
occur  in  a  perfect  fluid,  since  (5l)  would  imply  ^&   -   0   and 
then  all  quantities  '.^j-otTld  bo  continuous.   (IJ-.  cannot  vanish  by 
definition  of  a  phase  change  front.) 

Ill,   Geoinetrical  __Acoust i c s  -  Weak  Shocks 

A.   The  Variational  or  Acoustic  Equations 

Let  us  consider  a  one  para:Tieter  set  of  solutions  of  the 
equations  of  continuum  mechanics  dcpoiadinn;  upon  a  pararaetcr  h 
Wo  assunie  the  solutions  are  diff erentiablo  x-zith  respect  to  yi 
at  yi=   0   and  we  differentiate  (6),  (7),  ($^-(|?^>,  (i|.2)  with 
respect  to  yi  at  y^^  =  0,  obtaining 


(1)       p^  +  (piu  +  pu^)y^      =    0 

"i 


1  ^  .      0 


12)      u.      +U.U,      +   u.u.       -   -  T.  .      +   --4  T.  .        =   f. 

1.11  11  oil  d        1',  X 

t  ^       X  .  ^       X..       '  ''x  .       p  '^X  . 

«J  J  J  J 

/        N  •  •  •  1        r 

(3)       e .   +  u.e       +  \i,  e       +  --     q ,       -    T.  .u.        -     T-  ,u,.       ] 

^ -^  t  IX.  1    X.        p     ■   "1  111  11     1 

1  1      '  X .  ^      X .  -      X 


1  J  J 

D 

■    ^^i 

X  . 
1  J 


-   "^    [q.       -    T.  .u,       ]    =    0 
2    '■^1  ij    1 

0  X.  ■^      X. 


ik)      T.  .    =   T  .. 
^  ■  iJ  Ji 

(5)      e   =   Spp   +   e^s 

(d)      p  =  p   0   +  p„s 
P  2 

(7)     e  =  e^p  +  ej 

P  o 

In  these  equations  a  dot  denotes  differentiation  with 

respect  to  >^  .   Those  equations  are  linear  in  the  derivatives 
•   •   • 

p,  s,  u. ,  etc.,  which  arc  called  the  variations  of  the  corres- 
ponding undifferentiated  quantities,  or  the  acoustic  density, 
entropy,  velocity,  etc.   The  equations  are  called  the 
variational  or  acoustic  equations,  and  the  undifferentiated 
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quantities  arc  called  the  basic  flow.   It  is  to  be  noted  that 
a  solution  of  the  homogeneous  variational  equations  (with 
f.  =  0)  is  orivcn  by  the  x.  derivative  of  a  solution  of  the 
oriF;inal  equations  (with  f .  =  0)  . 

In  order  to  differentiate  the  discontinuity  conditions, 
vje  must  take  account  of  the  fact  that  the  discontinuity  surface 
may  depend  upon  r\  .   Thus,  from  II-(8)  we  have 


(S)  ^  +   f,  X. 


Since  only  normal  displace nents  of  the  surface  need  be 
considered,  vie   may  satisfj'-  (v6)  by  setting 

1 


(9)  x^  =  -^i'l'^^x.^' 


Making  use    of    (9),    we    obtain  from   (2©)-(aq) 

*       0       _T  /? 


(11)       [p^u.Au.  -  v.AT,.]-~[p,U,Z^u,-   v^AT,,]      vj(^2    ,-1/2  ^ 


(12)       [o^TJ-LA(e  + J  u^)    -  A(T.  .u^-  q.)v    ] 


-    fPlUi^(e-^  I  u2)  -  A(T.    u.  -  q    )v    ]^^^  vji^l  )-V2^ 


J 


K  k 


The  inequalities  j'-ield  no  definite  results  i\?hen  differentiated 
ujiless  the  equalities  hold  for  Y?=  0.   Then  thoj?-  also  hold  for 
the  derivatives,  unless  only  positive  yi    •'^^'^  considered,  in 
which  case  the  inequalities  '.lOld  for  the  derivatives.   Equations 
(10) -(12)  arc  linear  in  the  acoustic  quantities,  and  may  be 
called  the  acoustic  discontinuity  conditions. 

B.   Acoustic  Discontinuities 

'•men  the  basic  flow  is  continuous  (i.e.  p,  =  p.^j    '^^2_   ~  '^2' 
etc.),  (10)-(12)  simplify  to 

(13)   pv^Au^  -  -UAp 


11. 

Hk)      pUAu.    -  V  .AT.  . 

J  ""J 

(15)  pFA(e  +  u^.u^)    =   V  .  A(Tj[  .u^+ T^  .u^  -  q.)       . 

By  using   (111),    wg  may   simplify   (15)    to 

(16)  pUAe   =   V  .(T.  .  Au.    -  Aq  . )      . 

If  Tj  =  0  WO  have  an  acoustic  contact  discontinuity,  and 
(13) ,  [Ik) ,    (16)  become 

(17)  v^Au^  =  0 

(18)  V  .AT.  .  =  0 

(19)  V .  Aq.  =  0   . 

If  U  7^  0  but  ^q   -   0  we   have   an  acoustic   phase    change 
front  and    (13)    becomes    (1?)    v.rhile    (I4),    (I6)    become 

(20)  Zi>u^   =   (d-t)~^v.AT.  . 

(21)  TJ   =    (oAe)'^v^(T^^AA^    -   AQj)       . 

If  U  j^   0,  Z^  p  ^  0  wo  have  an  acoustic  shock  and  (Ik) 
becomes  ("lO)  »  which,  with  (13)  yields 


(22)  TJ  =  +  ._^J._....U   . 

''    A  p 

In  a  perfect  fluid  the  following  further  simplifications 
result: 

(23)  Z!s.p  =  0  contact  discontinuity 
( 2'-!-)      -^"i  ~  -(p"^')'  v.A.p      shock  or  nhase  front 


25)         LT  =  i  (  ^  shock 


c* 


12. 

(26)  IT  =  -(pAe)~  v.Aq-      phase  front 

(27)  pU9As  =  -v.Aq..  any  discontinuity 

The  last  equation  follows  from  II-(7),  II- ( 1+2  ^f'S'  i^  (l6). 

If  we  nov;  consider  a  perfect  fluid  without  heat  flo^^r 
(q.  =  G)  then  phase  fronts  become  impossible,  and  at  a  shock 

(27)  yields 

(28)  As  =  0   . 

TJslnp-    (28)    in   (5)    yields  Ap   -  P  Ap,    and   thus   frora   (25)    we 
have    for    the   velocity   of   an   acoustic    shock 

(29)  n  =  ±^-p-  =  i   c      . 


The  quantity  c  =  Jp      is  called  the  sound  speed,  and  it  may  vary 
with  position  and  time  since  it  depends  upon  the  basic  flow. 

If  the  definition  of  U,  (31),  and  II-(9),,  II-(IO)  are 
used  in  (29)  we  find 

(30)  (tt  +  ^^i'i^x  )^  =  ""^^^l    '      • 

This  is  a  first  order  partial  differential  equation  which  must 
be  satisfied  by  an  acoustic  shock  in  a  perfect ^  non-conducting 
fluid.   The  coefficients  depend  upon  the  basic  flow. 

If  (|),  =^  0,  we  may  write  the  equation  cj)  =  0  in  the  form 

(31)  t  =  W(x^)   . 
Then  (30)  becomes  an  equation  for  W 

(32)  (u  W^  -  1)2  =   c^l         . 

'^   J  '  J 

If  u.    =    0,    (32)    becomes    the    oiconal    equation 

(33)  \  =   ('v^    )^  =  W^    +i/    +\'jI        . 

c  J  12  3 
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If  u .  =  constant,  (32)  can  be  reduced  to  (33)  by  making  use  of 

•J  

the  change  of  variable  x.  =  x.-u.t.   The  equations  (32),  (33) 
a-PPly  only  when  c  and  u.  are  independent  of  t,  which  v;e  will 
henceforth  assuiTie  to  be  the  case. 

G.   Wave fronts  and  Rays 

The  acoustic  shock  surfaces  W  =  constant  v;ill  nox^j   be 
called  wavefronts  because  the  inte.cration  of  equations  such  as 
(32)  has  been  studied  in  georaetrical  optics  wnere  that  term  is 
used.   There  is,  of  course,  a  complete  theory  of  such  equations 
and  we  will  merely  specialize  it  to  this  case. 

First  we  define  the  wave  normal  p.  by 

{3h-)  Pi  =  w^^   . 

Prom   (32),    {3k-),   P-    satisfies    the   qiiadratic   equation 

(35)  ^^l^i  "  -"-^      ^  °~Pi      • 

At  each  point  in  space  (35)  defines  a  surface  containinfr.  the 
endpoints  of  the  possible  wave  normal  d. .   This  is  c^lDed  the 
surface  of  wave  normals.   The  reciprocal  of  the  radius  from 
the  oric-in  to  this  surface  in  any  direction  is  the  velocity  v 
in  that  direction  of  a  possible  vxavefront,  since  by  I-(IO), 
I-(31),  I-(32) 

(36)  V  =  (p^)-^/2   ^ 

By  examininp;  (35)  one   finds  that  the  wave  normal  surface 
is  a  siirface  of  revolution  about  the  direction  u.  .   Its  section 
by  a  plane  containing  this  direction  is  a  hyperbola,  parabola 

or  ellipse  according  as  the  basic  flow  is  supersonic,  sonic  or 

2  1/'' 
subsonic  (i.e.  accordinp-;  a?  u  =  (u.)  -^  ^  c).   If  u.  -  0  the 

-"■        _1      -'- 
wave  normal  surface  is  a  sphere  of  radius  c   .   If  Q   denotes 

the  anr-le  between  p.  and  u.  then  in  the  subsonic  case  there  is 

one  wave  normal  for  each  value  of  6,  and  the  wavefront  velocity 

v  =  c+u  cos  0.   The  same  applies  to  the  sonic  case,  except 

that  there  is  no  wave  normal  for  &  =  n. 


Ik. 

I.i  the  supersonic  case  there  are  two  solutions  p.  for 

-1  C  4- 

each  direction  if  C  <  0  <  cos   (  — )  ,  and  then  v  =  u  cos  *-)  1  c, 
the  upper  sign  applying  to  that  branch  of  the  surface  nearer 
the  origin.   V/hen  cos"  (^}  <  ©  <  cos"  (-  ^)    there  is  one 
solution  p,  for  each  direction  and  v  =  c  +u  cos  9.   If 
cos"  (— )  <  9  there  is  no  solution.   These  results  are  readily 
understood  in  terms  of  vectorial  addition  of  ix.  and  v.. 
The  Hamiltonian  function  H{x^,Pj^)  is  defined  by 

(37)  H(x.,p  )  =  c^'p.  -  (l-u..pj 

In  ter;ns  of  this,  certain  curves  x.((r)  ,  which  we  call  rays, 
are  introduced  by  the  equations 

dx  —^ 

(33)  ^  =  X^Ip,  =  ^\   (c^p.  +u.cjp^)   . 

?Iere  X  is  an  arbitrary  non-zero  factor.   From  the  last 

expression  in  (33)  we  see  that  the  rays  are  orthogonal 

trajectories  of  the  ;>ravc-c'ronts  if  and  only  if  u.  is  zero  or 

dx.  ^ 

■oarallel  to  p.,  since  then  -r-— ■  is  uarallul  to  p.  which  is  a 
- 1'  dcr    -  -^1 

wave  normal . 

Prom  the  general  theory  of  first  order  equations  it 

follows  that  the  rays  satisfy  Hamilton's  equations,  Lagrange's 

equations,  Fermat's  principle,  etc.   Furthermore,  by  means  of 

the  rays  vjc  can  construct  the  solution  of  (32)  corresponding 

to  any  given  initial  data.   The  v/avefronts  can  -also  be  found 

dx. 
by  Huyghen'  s  construction.   The  vector  -r;—  also  lies  on 

Fresnel' s  ray  surface,  which  is  reciprocal  to  the  wave  normal 

sLirface . 

D.   Variation  of  Shock  Strength  Along  a  Ray 

To  determine  the  variation  of  shock  strength  (i.e.  Ap? 
Ap,  etc.)  along  a  ray,  we  must  consider  the  acoustic  equations 
for  a  perfect,  non-conducting  fluid.   3y  specializing  (l)-(3)j 
(6)  or  by  differentiating  I-(I|-3)  ,  I-(;4.i;)  we  obtain 


,  ^c  ,;  i)    .    v">^  1  ;■ 


T'  .-' 


''■[•:-  i  ■    ■■; 
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(39)      p^  +  (pu^+  pu^) 


X. 

1 


0 


(JiO)       u.    +   u   U.       +    U.U.       +    D~    D^    -    po"    p.,      =    f. 


•  • 


(h-l)s,   +U.S+U.S  =0 

t  J    X.  J     X. 

ik?-)  P   =  PsS+  PpP      . 

We   now   consider    the    quc'.ntltios   u.  (x  .  ,V./(x  . )  )  =  u.  (x  . )  , 
p(x./-J(x  .)  )  =  p(x.)  ,    s(x  ./-/(x    )  )  ='g(x.)    and  p(x  .,>/(:,:.))=  p(x. ) 

JJ  JJJ  J  JJ  J 


•  • 


I'liese  arc  the  values  of  u.,  p,  s  and  p  on  a  wavGx'ront,  and 
there  are  two  such  functions  vjhich  we  define  in  this  way,  one 
corresponding  to  each  side  of  the  wavefront.   The  derivatives 
of  cither  of  these  functions  are  given  by 

u.    =  ti.   +  u.  W 
1       1      1 ,  X  . 
X.      ^^j     t   J 

^^3)  ^.  =  P^.  +  Pt  '^. 

J     J       J 

^_         •         • 

J    J      J 

Combining  (39)-(l-!-2)  with  ([|3)  yields 

(I1J4.)   p,[l-u.N   ]  +  u.  [-pVJ,^  ]  =  -pil.   -  p   u.  -  pu    -  p   u 
1      t      1        X.     1       X.     1 

'      1      t     -^   J  X 

"  ^-  ^>p"^px.-  Psp"\.^  pp'^'^x.-  X  ""r'-A 

^     1         1        ix.^'-^x. 

J       J 

-  p'^^Po   P  -  p'"^Po   "^ 
^x.       Px. 
1         1 
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One  such  set  of  equations  Is  obtained  on  each  side  of 
the  wavefront.   We  now  subtract  each  equation  on  side  tv;o  from 
the  corresponding  equation  on  side  one.   Then  remembering  that 
the  basic  flov;  was  assioi'ded  to  be  continuous,  and  that 
As  =  As"  =  0,  we  obtain 

(1+7)   [1-u.lC  ]  A(Pi.)  +  [-pV/„  ]A(u.  ) 

=  -pAu.   -  P.,  All.  -u.  ^p-u.A?^^^ 

X.      1  X.  1 

^      i  "^   j       t  i 

_2      -1       —     -1  —  _       _ 

=  [p  Px  "  P   ^o   ^Ap-P  P   Apx  -  ^i   ^u.-u.Au. 
i      Px ,       P       i    ^x .    -J    -^    X 

()49)   [1-u.W.^  ]A(s.)  =  -s„  /^u. 


J  J 


Equations  (!4.7)-('!-9}  are  a  set  of  five  linear  inhomogeneou£ 
equations  for  the  five  quantities  -Ap ,  ,  Au.   and  Z^  s^, .   Let  us 

G 

compare  these  equations  with  the  shock  conditions  which  are  a 
set  of  five  homogeneous  equations  for  Ap,  Au.  and  As. 
These  latter  equations  can  be  rewritten  in  the  form 


(50)  UAp   +   ov^Au^  =   0 


(51)  v^p  /^p   +   pUAu^   -I-   v^p^Ziks   =   0 

(52)  pUe„As   =   0 


1-^ 


s 


If  (1+7)  is  multiplied  by  -^  ,  (1+G)  by  -2-  and  (Il9)  by 


]\-'^  Jvt 


pe,,  ^"k  V  'k 

- — ~   then  the  coefficient  matrix  on  the  left  side  becomes 

identical  with  that  of  (50)-(52).   Therefore  ( [j.7 )-( 1-1-9 )  can  have 
a  solution  if  and  only  if  the  right  sides  of  these  equations 
satisfy  a  certain  orthogonality  condition,  since  by  assumption 
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the  homogeneous  equations  have  a  non-zero  solution.   To  obtain 
this  condition  we  eliminate  As^  from  (i|-8)  by  means  of  [k9)  > 
and  similarly  eliminate  As  from  (51)  t)y  means  of  (52).   The 
resulting  sets  of  four  equations  in  four  unknovms  vfill  both 
have  the  same  coefficient  ruatrix,  which  becomes  symmetric  if 
we  multiply  (i|7)  and  (50)  by  p   and  (I4.8),  (51)  hy  p.   Moiij  the 
solvability  condition  for  (l-!-7)j  (50)  is  that  the  rinht  side 
be  orthogonal  to  the  solution  of  the  transposed  homogeneous 
equations,  which  are  just  (SO)  ,    (51)  because  of  the  syiTuictrj. 
Thus  we  obtain  the  condition 

( 53  ;   -P^  (  P  Au,   +  P„  ^ u.  +  u .   Ap  +  u .  A  p,^  )  A  p 

1  1 


2,  r  -2       -1     -,  ^-      -1  ^  - 

+  P  (Ip   Pv  "  P   P^   jAp-p^p   Ap, 


^i      Px_.       P 

.-1 


1 


-  u.   Au.-  u.Au.   )Au.  +  U  po^s  Au.v.Au.  =  G 

^x.    ^        ^         -X.    ^      '  ^  ^i    ^  ^   ^ 
3  J  "^ 

Equation  (53)  can  be  considerably  simp].ified  a;:id  solved 
explicitly.   First,  v/e  note  that  A  u.=  Au.  and  Ap  =  Ap. 
Then  we  eliminate  A  u.  in  terms  of  -A  p  by  {2l\.)  ,    and  we 

eliminate  A  p  by  the  relation  Ap  =  c  Ap,  which  follows  from 
(5),  (23).   In  this  way  we  obtain  from  (53)  an  equation  in 
vjhich  the  only  acoustic  quantity  a"opearing  is  Ap.   This 
equation  is,,  after  dividing  by  ^p, 

.   u.(Ap)     2u  Ape,. 

v^-t-/     X^  IT  'x.     1      J      1       2         <i  -5 

IX.  X  .    c        C  C" 

1         1 

« 

p  v.Ap   v.(Ap)      X   1  J 
^x .  1   -     1    ■  X .     :: . 

1       +  -..       =■        -1 


>u3        " 


.  v.u.v.    ^   ^^   (4^)   ^p  s   V.  -^=  0 
1  J  i_.^_   -^2    u   ^plJ  'x^.   ^s  x^.  J  p^p 

0 
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Although  (5-0  appears  to  be  a  partial  differential 


equation  for  Ap,  it  can  be  reviritten  as  an  ordinary  differential 

eouation  alone  a  ray.   To  accomnlish  this  we  choose  \= 

in  (33)  so  that  the  equation  of  a  rajr  becomes  2c  7p  • 

dx.        u. 

(55)  ^  =  V,  -  -i   . 


dc 

u. 


Then  for  any  function  g  we  have  p,,  (v. )  =  4§.,  where  — §. 

"y^*     1     C        O-Vj  0." 

denotes  differentiation  with  respect  to  o"  along  a  ray.   Making 

use  of  this  fact,  and  observing  that  the  ninth  term  in  (pU) 

2 
vanishes  since  v.  =  1  and  thus  v.  v.    =0,  iSk)    becoraes 

X  . 
J 

(5°)   T7>  lO^  ^^ ■-^—  +  (  V. }   .  -  U.  (  — ) U.     =  O 

^-^     dCT   ^'    pc      1   ex.    i^c  X.    c    1 

J 

In  obtaining  (56)  from  iSk)    It  was  necessary  to  divide  by  — rr- 
and  to  make  use  of  (29)  in  the  form  11=  c. 

Equation  (56)  can  be  solved  at  once,  with  the  result 


(57) 


pc 


..2l 


or 


^o 


,  1  ,      11 

u .  ( — )  + u . 

1^  C   X.      C     1 

1  X.J 


do- 


er 

The  first  exponential  factor  exv      f  [-(v.)   ]d<r  has  a  simple 

'  4      ^  '^ 

geometrical  significance  xifhen  u.  =  0.   In  this  case  or  denotes 
arc  length  along  a  ray,  as  is  shoi^rn  by  {SB)-      Then  by  following 
the  method  of  Luneberg  (op.  cit.  ^o,  86)  we  find  that 

exp   j  (-v.)^  d<X  ~  T-   where  K  is  the  expansion  ratio  along  the 


<^o 


K 


ray.  This  expansion  ratio  is  the  limit  of  the  ratio  of  the 
normal  cross  sectional  area  of  a  tube  of  rays  at  (f  divided 
by  the  corresponding  area  at  <T  ,  as  the  area  tends  to  zero. 


J  >     1 
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Of   course,    the    tube    of   ra;-s    is    si^pposed   to    contain   the   ray   in 
question,    and   to   converge   to   it.      Thus   vrhen  \x^   =   G,    (57) 
becomes 


fr,q>  (Ap)"^    _   1 


LAPif 


o 


Since  ^=~^   is  the  normal  ivimr-  in  velocity  across  an  acoiistic 
pc 

shock  and  Ap  Au  represents  the  poucr  per  unit  ai'^ea  In  tne 
shock,  (53)  asserts  the  conservation  of  shock  front  energy  in 
a  tube  of  rays . 

Equation  (58)  yields  the  shock  strength  Ap  at  any  point 
on  a  ray  in  terms  of  its  value  at  one  point  on  the  ray.   Once 
Ap  is  known,  Ac  and  Au,  can  be  im^ricdiately  computed.   Since 
(53)  just  involves  K,  Ap  can  be  found  once  the  ra3'-s  are 
determined.   It  is  usefuJ.  to  note  that  the  normal  cross-section 
of  a  tube  of  rays  is  a  portion  of  a  wavefront,  and  that  the 
area  of  such  a  portion  is  proportional  to  the  Gaussian 
curvature  of  the  surface.   Thvis  K  is  al30  the  ratio  of  the 
G-aussian  curvature  at  a  ■'">olnt  on  a  wavefront,  to  the  corres- 
ponding curvatiare  of  another  wavefront  of  the  same  family  at  a 
point  on  the  same  vaj . 

E.   Reflection  and  Transmission  of  an  Acoustic  Shock 

at  a  Contact  Discontinuity 

I'lc   will  now   consider  the  reflection  and  transmission 
of  an  acoustic  shock  at  a  contact  discontinuity  in  the  basic 
flow.   An  incident  wavefront  v/ill,  in  general,  give  rise  to  one 
reflected  and  one  transmitted  v/avcfront,  and  to  a  modified 
motion  of  the  contact  discontinuity.   If  these  moving  surfaces 
are  represented  in  four  dimensional  x. ,  t  space  they  become 
stationary  three  dimensional  hyper surfaces .   They  are  joined 
together  on  the  contact  discontinuity  surface  \~^  ,    which  we 
will  represent  paramo trlcally  assuming  P  to  be  stationary 
by  the  equations 


(59)  x^  =  g^(A,,  w  ) 
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The    t>ireG  wave  fronts   ere    given  by 

(60)  w^  -  t   =   0      ,      I'J^  -  t   =   0      ,      i/  -  t   =   0      . 

Here,  and  below,  I,  R,  T  villi   denote  Incident,  reflected  and 
transmitted  respectively.   In  order  that  those  surfaces  be 
joined  together  on    r'   it  is  necessary  that  W  -Vf^M"^  on  P, 
that  is 

(61)  W^(g^)  =--   W^(g,)  =  -'^(g.)   . 

This  is  an  identity  in  ^,  v^  ,  and  may  therefore  be  dif forentiatod 
with  respect  to  ^,  yj  yielding 

(62)  (Wj   -  VTC  )g.   =  0   ,   (w3   -'^)?.„=0 

-r         m  -         m 

These  equations  show  that  the  difference  between  the  normals  to 

I      R  T       T 

W  and  VJ  ,  and  the  difference  between  the  normals  to  'J"   and  W 

are  parallel  to  the  normal  to  f^ . 

Ho  ^^fill  rewrite  these  conditions  after  introducing  the 

wave  normals  p  ,  p ',  and  p"",  and  the  vmit  norm.al  IT  to  P. 

Then  the  preceding  equations  imply 

■R     T 
p^  =  p-  +  a-jN 

(63) 

T     I 
p  =  p   +  a2N 

T    R    T 
Here  a-,,  a^   are  constants.   Prom  (63)  v/c  see  that  p  ,  o  ,  p 

and  F  all  lie  in  one  plane,  called  the  plane  of  incidence. 

Forming  the  cross  product  of  (63)  with  N  we  obtain 

(6L1)  p''  X  '^  =  P"  X  ?^  =  P^  X  '^   • 

If  wc  now  consider  magnitudes  in  the  above  equations,  and 
recall  that  the  length  of  a  wave  normal  is  the  reciprocal  of 
the  wavcfront  velocity  along  the  normal,  wc  obtain 
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(65)       —   sin  ©^  =  — -f-  sin  ©^  =  -i  sin  0^   . 
v^  v^  v^ 

Tho  first  eauality  is  the  la^^^  of  reflection,  and  yields  fe'^='it-© 
if  v'  =v"  5  wiiich  is  the  case  when  u.  =  0.   The  second  equality 

^  T 

is  Snell's  law  of  refraction.   The  angles  of  incidence  9", 

■R  .  T  I   ■p    T 

reflection  9"     and  transmission  G  are  octween  p  ,  p"',  p   and  IT. 

T 
T       -1  V        I  '"■' 

Since  9  =  sin  {—f   sin  P  ),  there  is  no  value  of  ©",  and  thus 

v^      I 
no  transmitted  wavcfront,  if  — =•  sin  ©  >  1,  which  is  called  the 

case  of  total  reflection. 

The  wavcfronts  and   P  divide  space  into  five  rcr^-ions  in 

the  neighbox-'hood  of  the  point  set  J  on  which  those  surfaces 

are  joined  together,   ";Jc  assume  that  as  J  is  approached  in  each 

of  these  regions  all  of  the  acoustic  quantities  have  finite 

limits  independent  of  the  direction  of  approach.   If  x-;e 

designate  those  rerions  as  ",1,2,3,14-  (see  figiire  1)  then  the 

limit  of  each  acoustic  quantity  will  be  indicated  by  a  si^per- 

script  corresponding;  to  the  region  in  vjhich  the  limit  is 

approached.   We  next  introduce  tho  jumps 

(Db)   Au^  =  u7  -u^  , 


AuJ 

0 

=  ^^i 

1 
-^i 

> 

^uf 

1       2 

=   Uj  -  u^ 

} 

AuJ 

2 

-Uf 

> 

a4 

-4-4 

> 

etc . 

It  is  assumed  that  A  u . ,  Ap  ,  Ap  ,  As   arc  knovm,  since 
these  are  the  jumps  across  tho  incident  wavefront,   3y  applying 

the  jump  conditions  we  wish  to  determine  the  other  jumps, 

1    ^ 
especially  the  rcflecter"  and  transmitted  ones  Au.  ,  A -^  in 

terms  of  those. 


From  tho  shock  condition   (2l|)  at  an  acoustic  shock,  we 


nave 


(6?)   Au°  =  -v^(pc)-^Ap''  ,   Au];  =   -v^(pc)-^Ap^-  , 

Au|  =  -v^^p'c')"^-^?^   . 


r  -,'•■;/ r 


^"V  ..       ;  -fJ 


/•/■" 
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Here    p,    c   designate   quarit5.tles   on   the    "incident''    side,    and   p'  , 
c'    on  the   '"'transmitted"'    side    of   the    contact   discontinuity. 
Also,    from   the    definitions    of   the    jump,    vie   have 


(68) 
(69) 


Ap°  +  Zip^  +  ^p2  +  ^p3  +  Ap^  =  0 
Au?  +  Au^-   +  Au?   +  Au^   +  A  u^  =   0 


Since  pressure  is  contin'ious  across  a  contact  discontinuity, 
Ap  =  Ap^  =  0,  and  (63)  becomes 


(7c) 


^\p°  +   Z^pl  +  Ap3  =  0   . 


Now,  we  take  the  component  of  (69),  normal  to  P  and  recall 
that  the  velocity  is  continuous  across  a  contact  discontinuity. 
Then  we  have  N^A  u^  =  IT^Au^  =  0  and  (69)  becomes,  after  using 
(67) 

(71)   (pc)"^cos  G^Ap°  +  (pc)'^cos  e^^p^ 

+  (p'c'  )"-'-cos  e^Ap'  =  C   . 

Equations  (70)  and  (71)  are  two  equations  for  the 
reflected  pressure  iurap  A"o   and  the  transmitted  pressure  jumt) 
-Ap  .   The  solution  of  these  equations  is 


(72) 


Ap  =  ^p 


T  T 

COS  0     cos  © 


pc 


P  c 


COS  0"^     COS  & 


R 


0  c 


pc 


(73) 


COS  0    cos  0 


-Ap-^  =  Ap° 


pc 


pc 


rn  -1 

COS  0   _  cos__02 
"p  c     "    ~P^^ 


The  coefficients  of  Zk p   in  (72)  and  (73)  are  the  reflection 
and  transmission  coefficients,  respectively. 
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If  u.  =  0  on  the  incident  side  of  tlic  contact  discontinuity 

■R      T  T?  T 

then  V   =  v'"  and  (65)  shows  that  Q"  =  n  -  ©  ,  provided  II  is 

oriented  so  that  0  <  G '"  <  4  .   Then  cos  9"  -   -cos  9      and  (72), 
(73)  become 

(71+)      Ap^  =  tt|^-^^° 

3  _   2    »   o    ^      _  _     DC  COS  ^'^ 


(75)     -Ap-'  -  T+z  ^"' 


7.  = 


p'c'  cos  9-'- 


Equations  ( 7[i.)  and  (75)  are  i'^resnel's  formulae,  usually  deduced 
for  reflection  of  periodic  plane  vjaves  from  a  plane  interface. 

I^.   Expansion  I^atio  for-  Straight  Rays 

The  expansion  ratio  '1   was  introduced  in  (53)  and 
determines  the  variation  of  the  shock  strength  along  a  rajT". 
It  may  be  defined  as  the  ratio  of  the  Gaussian  CLirvature  of  a 
wavefront  at  one  point  on  a  ray  to  that  at  another  point  on 
the  same  ray.   In  a  uniform  basic  flow  x-^ith  constant  sound 
speed  the  rays  are  otrai-Jiht  lines  ortho,e:onal  to  the  wavefronts. 
In  this  case  the  vario-as  vravcfronts  of  a  fa. .illy  of  wavefronts 
are  parallel  surfaces,  according  to  the  definition  of  parallel 
surfaces  employed  in  differential  geometry.   If  R-,  ,  Ro  are  the 
principal  radii  of  curvatu^r'e  of  a  surface  at  a  point,  then 
R,+Z,  Rp+Z  are  the  corresponding  radii  of  a  parallel  surface 
at  a  point  at  distance  z  away  on  the  same  orthogonal  (i.e.  ray). 
Thus 

-1  _    ^^r'2 

Prom  (76)  wo  see  that  K~   becomes  infinite  at  tvjo  points 
Z  =  -R  ,  -Rp  called  conjugate  points  on,  each  ray.   The  locus 
of  these  points  for  a  given  fa^'iily  of  wavefronts  is  called  a 
caustic  surface,  which  u.-^.ually  has  two  branches.   The  branches 
touch,  if  on  some  ray,  R-,  =  Rp.   A  point  at  which  the  two 
branches  of  the  caustic  sur-face  touch  is  usuallj?"  called  a  focus. 
Since  K   becomes  infinite  on  a  caustic,  the  acoiistic  shock 
strength  also  becomes  infintc  there. 


/--: 
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The  expansion  rat5.o  has  also  been  computed  for  a  iiniform 
basic  flow  v/ith  a  curved  contact  discontinuity  separcrtlnp;  two 
regions  of  different  constant  sound  speed.   In  this  case  the 
expansion  ratio  for  reflected  or  refracted  (i.e.  transmitted) 
wavef rents  becomes  important.   For  an  incident  spherical  wave- 
front  these  expansion  ratios  were  computed  by  J.  3.  and 
H.  B.  Keller,  "Determination  of  Reflected  and  Transmitted 
Fields  by  Geometrical  Optics",  Jour,  of  Optical  Soc.  of  America, 
Jan\is.vj   1950.   For  an  arbitrary  incident  wavef  rent  the  result 
is  much  more  complicated,  and  is  given  by  Stanley  Preiser  and 
Joseph  B.  Keller  "Determination  of  Reflected  and  Transriiitted 
Fields  by  Geometrical  Optics  II",  Research  Report,  Washington 
Square  College  Math.  Research  Group,  tlew  York  University,  1950. 
Some  parts  of  those  results  have  also  been  given  hj   others. 

G.   Example:   The  Shock  Tube 

As  an  example  of  the  utility  of  the  preceding  theory,  wc 
will  appl;/  it  to  the  shock  tube  or  one  dimensional  explosion. 
The  results  are  also  useful  in  understanding  some  of  the 
phenomena  in  spherical  explosions. 

Vi/e  suppose  that  all  quantities  depend  upon  x  and  t  only, 
and  that  onlj''  the  x  component  of  velocity,  vjhich  we  call  u, 
differs  from  zero.   The  region  -a  <  x  <  C  is  supposed  to  be 
filled  with  a  perfect,  non-conducting  gas  at  high  pressure, 
and  bounded  at  x  =  -a  by  a  rigid  wall.   The  region  x  >  0  is 
supposed  to  be  filled  by  a  jDcrfect,  non-conducting  gas  at  lower 
pressure.   As  the  basic  flow  i-ze  take  the  state  of  rest:   u  =  0, 
p  =  const.,  p  =  p-|  or  pp,  s  =  s-,  or  s^  according  as  x  >  0  or 
X  <  C.  (see  figujr'e  2)  , 

We  now  assume  that  at  t  =  0  the  acoustic  qu.antities 
vanish  for  x  >  0  and  that  for  -a  <  x  <  C  they  have  the  constant 

•     •     •     * 

values  p,  o,  s,  u  ==  0.   Thus  at  t  =  0  the  surface  x  =  C  is  an 
acoustic  discontinuity.   By  Huyghen's  construction  wc  see  that, 
if  this  gives  rise  to  acoustic  shocks  at  later  times  they  will 
be  located  at  x  =  c-,  t  or  -Cpt  where  c-,  ,  Cp  arc  the  sound  sj-^eeds 
in  the  basic  flow  for  x  >  0  and  x  <  0  rosnoctivcly .   Fur  the  i-mor  e. 
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the  strengths  of  these  t;hoc-:s  will  be  constant  since  the 
expansion  ratio  for  a  plane  wavefront  is  unity. 

To  solve  the  problem  we  assiime  that  the  sol^^tion  is 
piecewise  constant,  the  discontinuities  occurring  across  the 
shocks  and  across  the  contact  discontinuity  at  x  =  0,   The 
acoustic  quantities  to  the  right  of  the  shock  x  =  c-,  t  and  to 
the  left  of  the  shock  x  =  -Cpt  have  their  initial  values,  while 
between  the  shocks  the  acoustic  pressure  and  velocity  will  be 
denoted  by  p ,  u.   These  quantities  are  continuous  across  the 
contact  discontinuity  x  =  0,  and  we  assume  that  they  are 
constant . 

Prom  the  shod:  conditions  applied  at  each  shock  we  find 

(77)  L:i  =  u  =  -^     . 

P2C2  P^c^ 

Thus 

(pc)n 
(78)  P  =  P  -r 


pcTj+TpcTp 

^^^\  ^  =  P   (-p clTT-roTJ^       • 

The  above  results  (78),  (79)  will  apply  until  the  left 
shock  reaches  the  rigid  vjall  at  x  =  -a,  when  a  reflected  shock 
is  introduced.   Since  u  =  0  at  a  rigid  wall,  the  velocity  jurap 
across  the  reflected  shock  is  the  negative  of  that  behind  the 
incident  one.   Therefore  the  reflected  pressu.re  jump  is  equal 
to  the  incident  pressure  jiomp.   This  result  also  follows  from 
the  reflection  formula  (7'-!-)  when  p'c'  is  infinite.   The  solution 
with  one  reflected  shock  applies  until  the  reflected  shock 
reaches  the  contact  discontinuity,  when  it  is  partially 
reflected  and  partially  transmitted  according  to  (7k),  (75). 
The  reflected  shock  will  aga;i.n  be  reflected  at  x  =  -a,  and 
this  Process  will  continue  indefinitely. 

The  strength  of  the  first  shod:  traveling  to  the  right 
of  X  =  0  is  given  by  (78)  since  the  acoustic  pressure  in  front 
of  it  is  zero.   We  will  denote  the  strength  of  this  shock  by 
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c,p,  wliore  c=  (pc),[(pc)-,  +  (pc)p]~  .   Then  the  strength  of  the 
first  shock  traveling  to  'she  left  is  (a-l)p,  v;hich  is  also  the 
strength  of  this  shock  after  reflection  from  x  =  -a.   The 
shock  transmitted  through  the  contact  discontinuity  at  x  =  0 
will  have  strength  T(a-l)p  and  the  reflected  shock  v/ill  have 
strength  R(c-l)p.   Here  ^   and  T  are  the  reflection  and  trans- 
mission  coefficients  given  by  (7U)j  (7'?)  with  9"  =  0,  0^  -  -n: . 
In  this  \-jay  we  see  that  the  strength  of  the  n-th  shock 
traveling  to  the  right  of  x  =  0  is  R  ~  T(G-l)p  for  n  >  1. 

The  pressure  p   behind  the  n-th  shock  traveling  to  the 
right  of  X  =  0  can  easily  be  computed  by  summing  all  the 
"orossurc  jumps  across  th'.:  preceding  shocks.   In  this  way  we 
find  P-,  =  a  and 


(80) 


V      = 
n 


a  +  T(a-i; 


1-R 


n-1 


1-R 


D  =  2R   p 


n  >  2 


The  velocity/-  u  behind  the  n-th  shock  traveling  to  the  right 
of  X  =  0  is  then  given  by 

(31)  u. 


h 


Pn 
'l°l 


The  contact  discontinuity,  at  x  =  0  in  tho  basic  flow, 


moves  a  distance  d 


n 


-'~~   u   between  the  n-th  and  (n+Dr-'.t 
Cp   n  ^    • 


shocks.   Thus  the  total  distance  D  which  it  moves  until  the 


arrival  of  the  (n+l)st  shock  is 
(82)      D. 


n 


2a 
c^ 


PiC 


1^1 


na  .   ^1  n-1  -  f"  R^ 


\ 


— ^  (1-R  ,' 


P2^2 


Prom  (82)  we  see  that  D   oscillates  if  R  <  0,  and  D   varic 
Monotonically  with  n  if  R  >  0. 
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